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We recently proposed an empirical approach for the Landau-Pomeranchuk-Migdal (LPM) effects 
in photon emission from the quark gluon plasma as a function of photon mass. This approach was 
based on Generalized Emission Functions (GEF) for photon emission, derived at a fixed temperature 
and strong coupling constant. In the present work, we have extended the LPM calculations for sev- 
eral temperatures and strong coupling strengths. The integral equations for (f (p±)) and (g(p±)) are 
solved by the iterations method for the variable set {po, 1o, Q 2 , T, a s }, considering bremsstrahlung 
and aws processes. We generalize the dynamical scaling variables, xt, %l, for bremsstrahlung and 
aws processes which are now functions of variables po, go, Q 2 , T, a s . The GEF introduced earlier, 
<?t, 9t-i 9li 9l, are a l so generalized for any temperatures and coupling strengths. From this, the 
imaginary part of the photon polarization tensor as a function of photon mass and energy can be 
calculated as a one dimensional integral over these GEF and parton distribution functions in the 
plasma. However, for phenomenological studies of experimental data, one needs a simple empirical 
formula without involving parton momentum integrations. Therefore, we present a phenomenolog- 
ical formula for imaginary photon polarization tensor as a function of {qo,Q 2 ,T,a s } that includes 
bremsstrahlung and aws mechanisms along with LPM effects. 
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In this work, we present a study of Landau- 
Pomeranchuk-Migdal effects [l|, (LPM) in virtual 
photon emission from thermalized quark gluon plasma 
(QGP). The LPM effects on real photon emission from 
QGP have been reported 0,0] and an empirical approach 
in For the case of virtual photon emission in QGP, 
the processes that contribute at aa s order Q and the 
higher order corrections Q and LPM effects @ were 
well studied. In hard thermal loops (HTL) 0] method 
these processes occur at the one loop, two lop and higher 
loop levels represented by ladder diagrams. In the pho- 
ton emission calculations, the quantity of interest is the 
the imaginary part of photon retarded polarization ten- 
sor (SLP^). The dilepton emission rates are estimated 
in terms of this SJIP 1 ^, Bose-Einstein factor and Q 2 as 
given by Eq[TJ The 311^ including LPM effects is de- 
termined in terms of a transverse function f(p_i_) and a 
longitudinal part g(p±), as given by Eq[2] Q. 



For the case of virtual photon emission having small 
virtuality, the transverse vector function f(pj_) is deter- 
mined by the AMY equation (Eq[3]) and the longitudi- 
nal function by AGMZ equation (EqHJ) The en- 
ergy transfer function 8E(p±,p ,q ,Q 2 ,T,a s ) is given 
in Eq(5] The tilde represents quantities in units of Debye 
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mass, for details see [lOj . 
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I. GENERALIZED EMISSION FUNCTIONS 
FOR PHOTON EMISSION 



In the present work, we solved these Eqs. (|3l4p by itera- 
tions method at a fixed photon energy of yo/T=50. Alter- 
natively, these equations can also be solved by variational 
approach [ll|. In the following calculations, we have 
used two flavors and three colors. Using the iterations 



2 



method, we obtained p±.f (pj_), <?(p_l) distributions for 
different Po,qo,Q 2 , plasma temperatures (T=1.0, 0.50, 
0.25GeV) and strong coupling constants (a s —0.3Q, 0.10, 
0.05). We integrate these p_L.f (pj_), <?(p_i_) distributions 
to derive 7 T '^ as defined in the Eqs l6l71 The superscripts 
b, a in these equations represent bremsstrahlung or aws 
processes depending on the po value used. The subscripts 
T, L represent contributions from transverse (f(pj_)) or 
longitudinal (ff(pj_)) parts. l!^ a L are the quantities re- 
quired for calculating imaginary part of polarization ten- 
sor (see Eqf2]). Therefore, in the following, we empiricize 
these 4'° L . 
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In the remaining part of this work, we adopt the formulae 
and results of [10| presented at fixed T=lGeV, a s = 0.30, 
by suitably redefining the quantities for all tempera- 
tures and strong coupling constants. In Eqs l8l9llUl we 
define four dimensionless variables. The factor a s /0.3 
in above equations is required to match the definitions 
in present work with those of [l(J. The variable x\ is 
the real photon dynamical variable For virtual pho- 
ton emission from QGP, we define two more variables, 
x t,l given in Eqs lllll2"l l!^ a L are in general functions of 
{poj <7cb Q 2 i T, a s } and when plotted versus any of these 
jPQi Qo, Q 2 i they do not exhibit any simple trends. Fol- 
lowing jl0| , we define the generalized emission functions 
(GEF) g b T a L in Eq[T3] The GEF are functions of only x t ,l 
variables. These GEF (g^ a L ) are obtained from corre- 
sponding IrpL values by multiplying with c^ a L coefficient 
functions given in Eqs fTHfTSl The variable x in Eqs [TM^3l 
is xt for transverse part and Xl for longitudinal parts. 
The quantities xt,l an d c^ a L in Eqs llUfTSl are found by 
search for dynamical variables hidden in the solutions of 



AMY and AGMZ equations. 
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Figure [T] shows the results for GEF for bremsstrahlung 
(FigOJa)). The calculations are for a fixed photon energy 
(<7o/T=50.) but include six different cases of tempera- 
tures and coupling strengths mentioned in figure labels. 
The solid curve in (a) is the empirical fit to this emission 
function, given by Eq |19l 1 . The required c T coefficient 
function is given in Eq[l4l It has been observed that for 
low Q 2 , i.e., X2 < 2.0, transverse part of aws process be- 
haves similar to the transverse bremsstrahlung function. 
Therefore, we transform the low Q 2 transverse part of 
aws process as given by Eq|16l The resulting emission 
function is shown in Fig. 1(b). The solid curve is given 
in Eq |21[ which is same as solid curve in Fig. 1(a). 
The emission function for high Q 2 values (x 2 > 2.0) 
for transverse part of aws process is shown in Figure [21 
The c T and the emission function are given in Eqs ll5l201 
Similarly, Figures ((3{a,b) show the longitudinal com- 
ponents of GEF for bremsstrahlung (Fig. (a)) and aws 
(Fig.(b)). The coeffiicient functions and GEF are given 
in Eqs ll7l22ll8l23l These transformation functions are 
very complex. 2 . 



This fit given in Eq |19l is an improvement over the result reported 

2 The Eqs ll8l23l are slightly different from the corresponding equa- 
tions presented in Hoi . 
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FIG. 2: The dimensionless emission function gx(x) versus 
dynamical variable xt- The symbols are as in figure Q] Six 
different temperature and coupling constant values considered 
are mentioned in figure labels. The required c T coefficient 
function given in Eq.??. The solid curve is an empirical fit 
given by Eq |19l 



FIG. 1: (a) The dimensionless emission function g T (x) versus 
dynamical variable xt defined in Ea llll Six cases of temper- 
ature and coupling constant values considered are mentioned 
in figure labels in different colored symbols. The symbols rep- 
resent the integrated values of p± distributions as a function 
of {po, 1o, Q 2 , T, a s } values. These are transformed by c T co- 
efficient function given in Eq |14l Essentially, various symbols 
merge and can not be distinguished. The solid curve is an 
empirical fit given by Eq |19l (b) The dimensionless emission 
function g T (x) versus dynamical variable xt for X2 < 2.0. 
The transformation coefficients c T and empirical fit are given 
by Eqs fT5l2TJl 



II. GEF AND PHOTON RETARDED 
POLARIZATION TENSOR 



In the previous section, we used the results from the 
iterations methods to obtain the I^ a L values by integrat- 
ing the distributions. We transformed these into GEF 
(g^ a L ) functions shown in Figs[T][31 We fitted these by 
empirical functions given in Eqs ll9[T2"3l Using the em- 
pirical functions, for any po, qo, Q 2 , T, a s values, we can 



generate the I^ a L {x) values, circumventing the need to 
solve the integral equations. Thus, we have empiricized 
the l!p a L values in terms of GEF. Hence, using GEF and 

the c^ a L , the imaginary part of photon retarded polariza- 
tion tensor is calculated, as in EqCUflO]. 
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Here, the superscript i denotes {b, a} depending on the 
value of the integration variable po- 3 - We have calcu- 
lated imaginary photon polarization tensor and dilepton 
emission rates using Eq |24l and made a detailed compar- 
ison with the results of (8j. For this comparison, we 



The factor T in (Tmf,) in the Eq. [24] is arising from the tilde 
transformation. This extra T cancels the -jp- factor coming from 
tilde transformation of /, g functions. This T was missing in [id . 
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FIG. 4: (a,b) Dilepton emission rates using GEF method 
shown in symbols and compared with results of [12] repre- 
sented by blue lines. All the details are mentioned in figure 
labels and text. 



FIG. 3: (a) The dimensionless emission function g\{x) ver- 
sus dynamical variable xl defined in Eq |12l The symbols 
represent the integrated values of p± distributions as a func- 
tion of {po, qo, Q 2 , T, a B } values. These are transformed by 
c b L coefficient function given in Eq |17l The solid curve is an 
empirical fit given by Eq |22l The temperature and coupling 
constant values are mentioned in figure labels in different col- 
ors, (b) The dimensionless emission function gl(x) versus 
dynamical variable xl- The transformation coefficients cf, 
and empirical fit are given by Eqs ll8l2Tfl 



photon g /T=20.0 and a s =0.05 (see (a) ). The GEF re- 
sults were generated using T=1.0GeV. Similarly in Fig- 
ure (b) we show rates for go/T= 0.50 and a s =0.30 (in 
fig.(b)). The GEF results were generated at T=0.25GeV. 
The agreement of GEF method with lines is seen to be 
very good, except at the highest values of £ . This devia- 
tion is caused because for the longitudinal partin Eq[24j 
we used photon momentum When this is replaced 

with photon energy as shown in Eql25| the agreement 

of our results with 12] is very good in the full range of 
Q/T. In the remaining part of this paper, we use only 
Eql23 



generated reference results using the program provided 
by F. Gelis [II]. The agreement of the GEF method of 
Eq|24l with the results of [l^] was observed to be very 
good. As an example, we show the dilepton emission 
rates in Figure |U Figure shows the GEF results in sym- 
bols compared with the results of [l2| (blue lines) at a 
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We will present more results in a different way by defin- 
ing reduced quantities. After obtaining the SJIP^ ver- 
sus Q 2 ,qo,T,a s by using EqJ35J we define the reduced 
polarization tensor as and reduced Q rec i as in Eqs l26l2"7l 
The reduced polarization tensors are calculated for dif- 
ferent photon energies, different coupling strengths and 
temperatures. We plotted these results in balck circles in 
Figs l5l6l versus Q re d- For comparison, results from [jjj 
are shown in red symbols. The agreement of these two 
symbols is seen to be very good from low to very high 
photon energies, qo/T ~ 0.05 — 50.0. 



III. PHENOMENOLOGY USING 
GENERALIZED EMISSION FUNCTIONS 

In this section, we obtain the phenomenological fits to 
virtual photon emission rates from QGP. From the Fig- 
ures [321 it is clear that the reduced quantities depend on 
only two variables, i.e., instead of {Q 2 , qo, T, a s }, we need 
only {Qred, Qo/T} to generate n re ^_as in Eqj28l This ob- 
servation was already reported in 



{ >.yTI re£ / — F(Q r ed, rp ) 



(28) 



In the limit of Q red -> 0, F{Q rech f ) -> F (f ). To study 
this further, we use Eg 1251 to generate imaginary part of 
polarization tensor for various values of Q 2 ,qo, T, a s . At 
first, we generate the /mil at a very low Q 2 , Q ~ 10 _4 go, 
for various values of qo/T, a s . Using the results, we con- 
struct Q r ed, ^n re( j. The results are shown in Figure [7] by 
symbols labeled GEF method. The results for different 
a s ,T merge into a single curve in Fig(7] We fitted this 
data by suitable functions as given in Eqs l31l along with 
their parameters. These are two different fits, one for 
qo/T <200 and the other for qo/T >100, with an overlap 
between 100-200. These are represented by solid curves 
and labeled F in figure. The QTUd below q /T < 0.020 
is approximately equal to this function Fq, as given by 
Eql2lH However, it should be noted that at ultra soft 
photon energies, the present formalism needs corrections 

For the case of finite Q 2 , we made empirical fits by 
choosing a function given in Eq |301 In this function, 
A,B,C parameters are function of qo/T and are deter- 
mined by fitting the Q 2 plots for various qo/T. These 
parameters values for various qo/T are tabulated and 
are shown in Figure [H It is very important to have an 
empirical formula to generate A,B,C values. Therefore, 
these A,C,B parameters were fitted by different func- 
tional forms as given in Eqs l33l36l38l The parameters 



• GEF V= 
-F«U/I) 



l l I I I 



J 



: o GEF 


q/T=30.0 j/ 


r c Aurencheefi, 


/ (b) 


■ - F <w» 


Mill 



o GEF 




FIG. 5: QIl r 



S - / — for var- 



e d plotted as a function of Q re d 

ious photon energies (qo/T) mentioned in figure. The imagi- 
nary polarization tensor includes all contributions from trans- 
verse components of bremsstrahlung, aws, and also from the 
corresponding longitudinal parts. The black circles represent 
the GEF method in Eq |25l The red circles represent the re- 
sults of fl2l ] . The solid lines in violet color represent the results 
using EqfHCTl 



are different for different qo/T regions. Therefore, de- 
pending on the requirement, one can select the relevant 
parameter set to generate A,B,C values. 4 Using these 



4 It should be noted that the functional forms having difference 
of power law for these fits demand high precision of their pa- 
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FIG. 6: SjII re( i plotted as a function of Q re d for low qo/T 
values mentioned in figure. The details are as in previous 
figure [5] 
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FIG. 7: The reduced imaginary part of polarization tensor 
defined in Eg 1271 versus qo/T. We have taken as Q ~ 10 _4 5o. 
The solid curves are fits given in Eg 1311 The symbols repre- 
sent the results from GEF method using Eg 1251 
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formulae we get A,B,C coefficients and we get Fo(x) from 
Eq|3Tl We use these in Eqj30]to generate ^n re( j. These 
phenomenological results are shown by solid curves in 



rameters. Therefore, one should not truncate these parameters, 
especially the power exponents given by p,pi,P2- 
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FIG. 8: The A,B,C parameters versus qo/T. The curves 
represent fits by suitable functional forms in Eqs l33l 1361 l38l 
discussed in text. These are useful to generate H re d using 
Eq |30l Apparently, at high qo/T, these parameters are con- 
stant, however this is very misleading. The present fits gen- 
erate quite well the small variations of these parameters over 
full region. Good quality A,B,C fits are required because, the 
Q 2 plots are sensitive to these parameters ane there is delicate 
cancellation of various terms in Eg 1301 



A(x) = a + {bx Pl - cx P2 ) for x < 0.1 

a = -4.84554727516 
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In conclusion, the photon emission rates from the 
quark gluon plasma have been studied as a function 
of photon mass, considering LPM effects at various 
temperatures and strong coupling strengths. We defined 
generalized dynamical variables xt,xl for transverse 
and longitudinal components of bremsstrahlung and aws 
mechanism. In addition, we defined generalized emission 
functions (GEF) namely <7t {-^t ) j9t ( x t ) {^l) :9l 
We have obatined empirical fits to these GEF. In terms 
of the GEF, we have calculated the imaginary part of 
retarded photon polarization tensor as a function of 
photon energy and mass, plasma temperature and strong 



coupling strengths. For phenomenological applications, 
we fitted the reduced imaginary polarization tensor by 
simple functions and provided necessary parameters. 
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